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We consider the critical behavior for a string theory near the Hagedorn temperature. 
We use the factorization of the worldsheet to isolate the Hagedorn divergences at all genera. 
We show that the Hagedorn divergences can be resummed by introducing double scaling 
limits, which smooth the divergences. The double scaling limits also allow one to extract 
the effective potential for the thermal scalar. For a string theory in an asymptotic anti-de 
Sitter (AdS) spacetime, the AdS/CFT correspondence implies that the critical Hagedorn 
behavior and the relation with the effective potential should also arise from the boundary 
Yang-Mills theory. We show that this is indeed the case. In particular we find that the free 
energy of a Yang-Mills theory contains "vortex" contributions at finite temperature. Yang- 
Mills Feynman diagrams with vortices can be identified with contributions from boundaries 
of moduli space on the string theory side. 
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1. Introduction 

Since the early days of string theory, it was observed that the free string spectrum has 
a density of states which grows exponentially with energy, and that the partition function 
Z = e~ l3H of a free string gas at a temperature T = 4 would diverge when T is greater than 
some critical value Th [00:11] • The Hagedorn divergence occurs for all known (super) string 
theories with spacetime dimensions greater than two. The physical meaning of the critical 
temperature Th and of the divergence has been a source of much discussion since then. 

In the late eighties, a few important observations were made which suggested that the 
Hagedorn divergence signals a phase transition, analogous to the deconfmement transition 
in QCD @,[5|||0] ■ At the Hagedorn temperature Th the lowest winding modes (with wind- 
ing ±1) around the periodic Euclidean time direction become marginal operators in the 
worldsheet conformal field theory Sathiapalan and Kogan |^]|| argued that above 

the Hagedorn temperature, the winding modes would condense in a fashion similar to the 
Kosterlitz-Thouless transition in the X-Y model and the worldsheet theory will flow to a 
new infrared fixed point. From the spacetime point of view, these winding modes (with 
winding ±1) correspond to a complex scalar field (ft living in one fewer spacetime dimen- 
sion (i.e. not including Euclidean time). Near the Hagedorn temperature, the spacetime 
effective potential for (ft can be written in a form 

V = ml((3)r<ft + A 4 <7 s 2 (</>» 2 + A 6 S s 4 (</>» 3 + • - • , m% ((3) oc T H - T . (1.1) 

If A4 is positive (negative), the phase transition would be second order (first order). In |7| 
Atick and Witten argued that for a string theory in asymptotic flat spacetime the transition 
should be first ordeiQ (i.e. A4 < 0) due to the coupling of the thermal scalar to the dilaton. 

While the one-loop Hagedorn divergence has been extensively discussed in the past (see 
e.g. Pjnj for reviews), Hagedorn divergences from higher genus amplitudes have been 



investigated rather little. In this paper we use a factorization argument to extract Hagedorn 
divergences for higher genus amplitudes. We show that they can be resumed by introducing 
various double scaling limits, which smooth the divergences. The double scaling limits also 
allow one to extract the effective potential ( |1 . 1| ) to arbitrary high orders. That a double 
scaling limits might exist for higher genus Hagedorn divergences was speculated earlier 
in 



rjJI and further discussed in |Tj| in a toy model motivated from AdS/CFT. 



1 That the transition is of first order can also be argued from the non-perturbative instability 
of the thermal flat spacetime discovered in ||. 



1 



Our discussion further highlights that Hagedorn divergences signal a breakdown of 
string perturbation theory due to appearance of massless modes and do not imply a limiting 
temperature for string theory (|,[5],|7| . 

The discussion of this paper will be rather general, e.g. applicable to string theories 
in asymptotic anti-de Sitter (AdS) spacetime. The AdS/CFT correspondence then implies 
that the critical Hagedorn behavior from high genera and the relation with the effective 
potential should also arise from Yang-Mills theories. We show that this is indeed the 
case. In particular we find that the free energy of Yang-Mills theory contains "vortex" 
contributions at finite temperature. Yang-Mills Feynman diagrams with vortices can be 
identified with contributions from the boundary of the moduli space on the string theory 
side. 

The plan of the paper is as follows. In section 2 we first review the one-loop result 
and discuss the physical set-up of our calculation. We then extract the critical Hagedorn 
behavior from higher genus amplitudes and show that one can find terms in ( |1 . 1| ) by 
defining suitable double scaling limits. In section 3 we turn to Yang-Mills theory. We 
discuss the structure of the large iV expansion for the partition function of a Yang-Mills 
theory at finite temperature and isolate the critical Hagedorn behavior. We conclude in 
section 4 with a discussion of some physical implications. 



2. High-loop Hagedorn divergences in perturbative string theory 

2.1. Review of one-loop divergence and set-up 

Consider a string theory consisting of a compact CFT times IR 1 '^. The one- loop 
free energy of the system at a finite temperature can be computed by the torus path 
integral with a target space in which the Euclidean time direction is compactified with 
period (3 = h and with anti-periodic boundary condition for spacetime fermions ||13|| . 



The Hagedorn singularity appears when the lowest modes with winding ±1 around the 
compactified time direction become massless ||[|]||. More explicitly, the mass square can 
be written as 

(^) 2 -(^) 2 ™ 

where the first term is the winding contribution and cq is the zero point energy of the 
string (in the winding sector). The second equality of ( |2~T| ) should be considered as a 
definition of the Hagedorn temperature. From (PTT|), m|(/3) — ► as (3 — ► (3h and becomes 



tachyonic when f3 < In spacetime, the winding ±1 modes correspond to a complex 
scalar field <p living in one fewer spacetime dimension (i.e. spatial part of the spacetime), 
which is often called the thermal scalar in the literature. We will follow this terminology 
below. We will also collectively call modes with general winding numbers (and no internal 
excitations) winding tachyons. Equation ( |2.1| ) applies to both bosonic and superstring 
theories with possibly different Co for different theories. 

The critical behavior of the one-loop free energy F\ as (3 —>■ (3h is controlled by that 
of the thermal scalar 

F x = -2 x i log(-V 2 + mj(/?)) + F finite , (3 -+ j3 H (2.2) 

where V 2 is the Laplacian on the spatial manifold. If the gap of V 2 along the compact 
CFT directions is bigger than rr£{f3), the singular part of (2^2) can be further written as 



{ml{(3))i d odd (2 ' 3) 

(m 2 ((3))z logrri^(P) d even 

F\ has a branch point singularity at m|(/3) — for all d. In particular, for d = it is 
logarithmically divergent as (3 — > /3jj 

Fx = -log(/3- p H ) + finite . (2.4) 

The above discussion should also apply to a static curved spacetime, for example, an 
AdS spacetime, even though an explicit computation of the one-loop free energy is often 
not possible. For an AdS spacetime, since the Laplacian has a mass gap, we expect the 
free energy for a thermal gas of AdS strings should behave as (|2.4|) when the Hagedorn 



temperature is approached (see e.g. |14| for further discussion). 

In this paper we will focus our discussion on d = or more generally those spacetimes 

(including AdS) in which ( |2.4| ) is satisfied, for the following reasons: 
1. The thermal ensemble cannot be defined in an uncompact asymptotically flat space- 
time due to Jeans instability. To make the canonical ensemble well defined, an Infrared 
(IR) cutoff is needed. One particularly convenient (and well-defined) IR regulator is 



to introduce a small negative cosmological constantB. For our discussion below the 
precise nature of such a regulator will not be important as far as it makes the thermal 
ensemble well defined. Such IR regulators introduce a gap in the Laplacian V 2 , which 
will be kept fixed in the limit T — > Th and thus will be greater than m 2 when the 
temperature is sufficiently close to Th- 
2. The Hagedorn singularity is sharpest at d = 0. While the free energy is singular at 
f3 = (3h for all dimensions, it is divergent only for d = 0. 

The logarithmic divergence of ( [2.4] ) at f3 — > Ph implies that the string perturbation 
theory breaks down before (3 = (3h is reached. Thus it is not sufficient to consider only the 
one-loop contribution to the free energy and higher genus contributions could be important. 
Below we will show that as (3 — > /?#, it is necessary to resum the string perturbation theory 
to all orders. We will then show that one can extract the spacetime effective action for the 
thermal scalar from the resumed series and that the divergences are smoothed out. 



When A4 in ( |1 . 1| ) is negative, i.e. when the transition is first order, there exists a 
lower temperature T c < Th, at which the thermal gas of strings becomes metastable. At 
a temperature T c < T < Th, the thermal gas is still perturbatively stable. Here we are 
interested in probing the critical behavior in the metastable phase (or superheated phase) 
as T — > Th from below. 

2.2. Higher loop divergences 

We now examine higher loop divergences as j3 — > (3h- For simplicity we will restrict 
our discussion to bosonic strings. We expect the conclusion to hold for superstring theories 
as well. 



2 In an asymptotic AdS spacetime it is possible to define a canonical ensemble in the presence 
of gravity, as discovered by Hawking and Page [15]. Hawking and Page also found that the system 



undergoes a first order phase transition at a temperature Thp from a thermal gas in AdS to 
a stable black hole. Treating an AdS spacetime with a small cosmological constant as an IR 
regularization of the flat spacetime, it is natural to identify the first order phase transition argued 
by |?J with the Hawking-Page transition. Note that the flat space limit, which corresponds to 
keeping g s small, but fixed and taking the curvature radius of AdS to infinity, is rather subtle. In 
this limit the stable black hole phase in AdS disappears and the Jeans instability should develop 
at a certain point. Also note that in the flat space limit, the Hawking-Page temperature goes 
to zero, which is consistent with the observation that a hot flat spacetime is non-perturbatively 
unstable at any nonzero temperature [pi. 
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The genus-g contribution F g to the free energy is obtained by integrating the single 
string partition function on a genus-g surface over the moduli space M. g of such surfaces. 
The potentially divergent contributions to F g arise from the integration near the boundary 
of the moduli space. 



Fig. 1: An example of a degenerate genus-6 Riemann surface. Each blob represents 
a surface of certain genus and thin lines connecting blobs represent pinched cycles. 




Fig. 2: Degenerate limits of a genus-2 Riemann surface. 

The boundary A g of Ai g is where a Riemann surface degenerates, which can be 
described by pinching cycles on the surface (for reviews see e.g. |!6| , [l7yiOH ) ■ There are two 
types of basic degenerations depending on whether the pinched cycle is homologous to zero 
or not. If the pinched cycle is homologous to zero, a surface of genus g degenerates into 
two surfaces of genus g\ and g<i (g = g± + g^) which are joined together at a point. If the 
pinched cycle is not homologous to zero, a genus g surface degenerates into a surface of 
genus g — 1 with two points glued together. One can pinch more than one cycle at the same 
time as far as they do not intersect with each other. On a genus g surface, the maximal 
number of nonintersecting closed geodesies is 3g — 3, so one can pinch at most 3g — 3 cycles 
at the same time. See fig. [I] and fig. |^ for examples of degenerate limits. 

Let us now examine the contribution to F g from boundaries of moduli space. The 
pinching of a Riemann surface can be described in terms of cutting open the path integral 
on the surface. The pinching is a local operation and so is cutting the path integral 

5 



(other than possible constraints from the zero mode integration). We follow the standard 
procedure as described in [ItJ. One has 



( 1 ) 9 = Yl<l hi t i (A(z 1 )) gi (Mz 2 )} g2 (2.5) 

i 

and 

i 

for the two types of basic degenerations, where (• • •) denote worldsheet correlation func- 
tions on a genus g surface and % sums over a complete set of intermediate states, q can 
be considered as the complex coordinate transverse to the boundary with q — > corre- 
sponding to the degeneration limit. Integration of ( |2.5| ) and (|2.6| ) near q — > yields the 
propagator 

G = Z} a /(_ V 2 + m ?) • (2 - 7) 

The contribution to the free energy from boundaries of moduli space can be extracted 
from diagrams like the ones in fig. [I] and fig. 0. One can treat blobs (representing surfaces 
of certain genus with some insertions) as effective vertices and thin lines (pinched cycles) 
as propagators. For (3 — > (3h and assuming that the spatial Laplacian operator —V 2 has a 
gap, then the propagator ( |2.7| ) for a pinched cycle is potentially dominated by that of the 
thermal scalar! , 

Since one can pinch at most 3<? — 3 cycles at the same time, naively we may conclude 
from ( p.8| ) that F g diverges as m-p^Ag-s for g > 2 as /3 ^ (3h- However, there are global 
constraints due to winding number conservation at each blob of fig. [TJand fig. [2| As a result, 



3 Note that it is not immediately obvious that the thermal scalar (or other winding modes 
along the Euclidean time direction) appears in the intermediate states from the point of view of 
calculating the free energy of a finite temperature string gas, since they do not correspond to 
spacetime physical states. Indeed in the one-loop calculation, they appear only after a modular 
transformation. However, it is clear that they should appear in the intermediate states from the 
point of view that we are working with a string theory compactified on a circle with anti-periodic 
boundary condition for fermions. 
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not all propagators can have the nearly-massless thermal scalar propagating through them. 
We will now show that the most divergent terms at genus g are proportional to 

' 9 > 2 • (2.9) 



Let us consider a generic degenerate limit of a genus g surface as shown for example 
in fig. |l|. Denote l/( n > m ) the number of vertices with genus n and m insertions. Then the 
total number L of pinched cycles (propagators) and the genus g of the whole surface can 
be written as 

2L = J2 mV(n ' m) i # = 1 + 5Z(y + n-l)v {n > m) . (2.10) 

n,m n,m 

The second equation of ( |2.10| ) can be obtained from the degenerate rules stated earlier. 
Alternatively, one can associate each insertion with a factor of g s and the total power of 
g s should be 2(g — 1). It is also convenient to introduce 

V = V {n ' m \ 9a = J2 nV(n,m \ (2-11) 

n,m n,m 

where V is the total number of vertices, g a is the apparent genus of the diagram (i.e. the 
sum of the genus of each vertex). Equations ( p. 10|) and (|2.11|) lead to 



L-(V-l)=g-g a . (2.12) 

Since winding numbers carried by propagators have to be conserved at each vertex, equa- 
tion ( |2.12| ) implies that the total number of independent windings in a diagram is g — g a . 
The maximal number of independent windings among different degenerate limits is then 
g, in which cases each vertex has the topology of a sphere. 
From ( |2.iq ) and ( |2.11| ), we also have 

V = \L-\Y,{m-A)V^ (2.13) 



and (|27LJ) and (|27L2]) lead to 



L = 2(g - 1) - ^(2n + | - 2)V^ . (2.14) 

n,m 

Equation ( |2. 14j) implies that the maximal number of propagators (pinched cycles) in a 
degenerate limit is indeed 3g — 3, obtained when only l/(°> 3 ) is nonzero. However, it is 
impossible to have all 3^ — 3 propagators to be divergent at the same time, i.e. to have 
all windings to be ±1, since by winding number conservation if the windings of two of the 
propagators coming out of a 3-point vertex are ±1, then the third one can only be 0, ±2. 
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Fig. 3: Two possible degenerate limits of a genus-3 Riemann surface which give rise 
to most divergent contributions. Each propagator has the thermal scalar running 
through it. 



Since at least one of the propagators going out of a 3-point vertex must have winding 
\w\ 7^ 1, if our purpose is to find the maximum number of propagators that can have 
w = ±1, one can ignore such a propagator. This implies we only need to consider those 
degenerate limits in which effective vertices have at least four insertions, i.e. m > 4. In 
the absence of V^°' 3 \ equation Q2.14 ) implies that 



L< 2(0-1) 



(2.15) 



where the equality holds when 



y(°> 4 ) 0, otherwise V {n > m) = 



(2.16) 



Thus we have proven that the most divergent term is of the form ( |2.9| ). See fig. |3| for 
degenerations which give rise to the most divergent contributions at genus 3. 

To summarize, the most divergent contributions at each genus have the following 
diagrammatic structure: 
1. Each vertex has the topology of a sphere and has four winding tachyon operator 

insertions with winding numbers 1,1,-1,-1 respectively. The total number of vertices 

in a genus g diagram is g— 1. The path integral over each vertex gives rise to an effective 

coupling 

g% = ( V+i(0)V+i(l)V_i(oc) / dhV-^z)) (2.17) 



Note that at (3 = Ph, the vertex operators V±i for the thermal scalar are marginal 
and ( |2.17| ) is well defined. Also A4 is (^-independent. 
2. The propagators are given by that of the winding tachyon ( |2.£| ). The total number of 
propagators is 2(g — 1). 
Thus the most divergent contribution to the free energy at genus- (n + 1) has the form 

(\ 2n 
^p)J k (^£f (2 - 18) 

where a n is a combinatoric numerical factor depending on the specific geometric structure 
of boundaries of moduli space. Determining these numerical factors from direct worldsheet 
computation is a rather complicated mathematical question, which goes beyond the scope 
of this paper. In the next subsection we will determine them using an indirect argument. 

2.3. Double scaling limits and the effective thermal scalar action 

In the last subsection we showed that the leading order Hagedorn divergences at all 
loop orders can be written as 



^9s , / fi'gAz 



7J 



with 



F s ing = - log(/3 - Ph) + «1 H 1" O n I ) + '" ( 2 " 19 ) 

mi \ m, ' 



Equation (|2.19|) suggests a double scaling limit 



P-Ph^0, 9s ^0, 13 ^ = finite (2.21) 

9s 

in which case all higher order terms in the series become equally important and we need 
to be resumed. 

How do we interpret the free energy F obtained by resuming the series? A clue comes 
from the structure of the degenerate diagrams summarized at the end of the last subsection, 
which resemble the Feynman diagrams of a |</>| 4 theory (see e.g. fig. |3|). Indeed the free 
energy of a |</>| 4 theory gives an asymptotic expansion which is precisely of the form ( |2.19j ) 
with specific values for the numerical coefficients a n . Given that string theory should 
reduce to a field theory in the low energy limit, and that here we are essentially isolating 
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an effective theory for the nearly-massless thermal scalar, it is natural to conjecture that 
( |2.19|) can be written as 



F sing =\og J dcf>d^e-<^-^ { 



2 \ 4 \ 2 (2.22) 

-log03-^)-2^ + 10^i + ... 



with is a c-number. Equation Q2.22|) determines a n to all orders uniquely and implies 



the following effective potential for the thermal scalar 

V = m 2 M*+W s (<P<P*) 2 + -- - • (2-23) 



In the next section we will show that the effective action ( j2.23|) and ( 2.22 ) arises from the 



critical behavior of Yang-Mills theories near the Hagedorn temperature. Using AdS/CFT 



this would serve as a proof of (|2.22| ) for string theories in an asymptotic AdS spacetime. 



Furthermore, since the factors a n in (|2.18|) and (|2.19| ) depend only on the mathematical 



structure of the moduli space of Riemann surfaces and not on the specific string theory, 
the Yang-Mills theory results serve as an indirect proof of ( |2.22j ). 

It is clear from equation ( |2.22| ) that Hagedorn divergences at each genus order in 



( |2.19| ) simply signal breakdown of the asymptotic expansion in g s due to that </> becomes 
massless. The — > limit is apparently smooth in the resumed integral expression 
fl2.22|) . When A4 is positive, i.e. when the transition is second order, the integral (|2.22|) is 



finite and non-perturbatively defined. For negative (or zero) A4, i.e. when the transition 
is first order, the integral ( |2.22| ) is not defined non-perturbatively and higher order terms 
in the effective potential are needed. In either cases the m§ — > limit is well-defined. 

Equation ( |2.22| ) implies that a n ~ nl for n large. This is in contrast with the (2n)! 
growth of the asymptotic behavior for the full free energy. Here we are only looking at 
contributions from boundaries of moduli space, which accounts for the slower growth. 
When A4 < 0, one can formally resum the series ( |2.19| ) or the second line of ( j2.22| ) using 



Borel resumation and one finds that the free energy contains an imaginary part of the 

_ 1 

form e due to the n\ growth of a n . Such an imaginary part can be interpreted as 
the tunnelling rate from the metastable thermal string gas to the true non-perturbative 
minimum (see also discussion in [|12])- 
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Here we have been focusing on the lowest spacetime modeu of the thermal scalar, 
which gives the most divergent contribution to the free energy. This explains the finite- 
dimensional integral in ( |2.22j ). From general covariance it seems natural to generalize 
( p.23|) to include derivatives 

S = J d d x^g (m 2 + ml<M>* + WA^f + ■■■) • (2-24) 

where d d x integrates over the spatial directions^. 

Let us now consider the generalization of the above double scaling argument to extract 
higher orders terms in ( |2.23|) . From equation (|2.10|) the leading contribution of a generic 



degenerate surface to the free energy can be written in the form 



2,-2 E„, fe ^"' 2fe) (2n+2 fc -2) 

U " - fl " (2.25) 



where in writing down Q2.25|) we have assumed that all propagators in a degenerate diagram 



carry winding numbers^ ±1 and that each vertex contains an even number of insertions 

m = 2k, k = 2, 3, ■ • -, due to winding number conservation. Now consider the double 
scaling limit 

(/3 ~f U) = finite, g s - (2.26) 

9s 

under which (|2.25|) is proportional to gjf with K given by 

oo oo 

k = j2Ys v(n,2k) ( 2n+2k - 2 - ka ^ ■ ( 2 - 27 ) 



n=0 k=2 



For a < 1, we always have K > for any choice of V^ n ' 2k \ At a = 1, we get K — for 
diagrams with y(°> 4 ) ^ only while K > for all other diagrams. In the double scaling 
limit ( |2.26| ) only the contributions of diagrams with K = survive. These are the most 



4 Recall that we assume that the Laplacian of the spacetime manifold has a mass gap. 

5 Note that for an AdS with a small negative cosmological constant, ( |2.24| ) applies to regions 
in the interior of the spacetime, since in AdS gtt component of the metric is nontrivial and the 
thermal scalar always has a large mass near the boundary. 

6 If there is a propagator carrying a winding number other than ±1, we can treat the two 
vertices connected by this propagator as a single effective vertex. Keeping doing this we obtain a 
degenerate diagram whose propagators only carry winding numbers ±1. 
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divergent contributions we isolated in ( p.!9|) and lead to the effective action (|2.23|) . Now 



let us set by hand A4 = 0, then in ( p. 27 ), l/(°' 4 ) = 0. The most divergent contributions 

4 

3' 



in the remaining diagrams are isolated by taking a = ~, at which K = for diagrams 



with V 0,6 7^ only and K > for all the rest. In other words now the most divergent 
contributions to the free energy can be written as 

F = — log(/3 - fa) + ° ig2s 3 + • • • + Cn ^ n 3 „ + • • • (2.28) 

which implies the effective potential 

V = mJ</></>* + A 6 (</></>* ) 3 + --- (2.29) 



where \q is related to the genus-0 six-point function of the vertex operators for the thermal 
scalar on the worldsheet. Now restoring A4 and combining ( |2.22j ) and ( 2.29 ) we would 
conclude that the effective potential can be written as 

V = ml<M>* + A 4 (#*) 2 + A 6 (#*) 3 + • • • (2.30) 

The same procedure can then be repeated to the next order by first setting A4 and A6 to 
zero and then extracting the most divergent terms in the remaining diagrams. One can 
continue this to arbitrary orders in (</></>* ) n and we find the effective potential^ 

00 

V = m%W + J2 + • • • (2-31) 

k=2 

The A2/C term is obtained by setting all vertices with m < 2k to zero and performing the 



scaling P-P H ~ gT l \ i.e. a = 2(1 - ±) in ( ggg) . 



Finally let us consider how to define various Ae, As, • • • from string amplitudes. Recall 
that A4 can be obtained from ( |2.17|) and ( |2.20|) . Naively one might want to define A2&; 
for k = 3,4, ■■• by the tree-level amplitudes of k winding 1 and k winding —1 modes. 
However, from factorization argument, these amplitudes are divergent at mi = 0. The 
divergences come from diagrams containing lower order vertices \2k' with k' < k and <p 
in the internal propagators, which can be found from standard Feynman diagrams for the 
action m\<j>4>* + Y^k~=2 ^k9s k ~ 2 ( ( P ( l ) *) k • ^2k is thus given by the sphere amplitude of k 
winding 1 and k winding —1 modes with the divergent parts subtracted. 



Note that the procedure is not well adapted to resum divergences due to vertices with genus 
n > 1. From ( |2.27| ), to have K = for n = 1, we need a = 2, in which case all genus 1 vertices 
with arbitrary number of insertions contribute equally. To have K = for n > 1, we need a > 2, 
then from ( |2.27| ), diagrams with large k become more dominant regardless of the value of n. 
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3. Hagedorn behavior from YM theories 

Our discussion in the last section was rather generic. In particular it should apply to 
type IIB string theory in AdS 5 x S5 or other string theories in asymptotic AdS spacetime. 
In an AdS spacetime with curvature radius R much bigger than the string and Planck 
lengths, there is a first order Hawking-Page transition at temperature Tjjp ~ much 
below the Hagedorn temperature Tjj ~ —k= at which the thermal string gas in AdS 

VOL 

becomes perturbatively unstable |y|. The discussion of the last section describes what 
happens if one stays in the superheated thermal AdS phase above the Hawking-Page 
temperature all the way to the Hagedorn temperature. From the critical behavior at the 
Hagedorn temperature one can then map out the potential for the thermal scalar. Aspects 
of the Hagedorn transition in AdS have been discussed in |18|,|14 . 



Hawking and Page's semi-classical discussion applies to IIB string theory in AdS with 
a cosmological constant small compared to the string scale and to the Planck scale, which 
corresponds to Af = 4 super- Yang-Mills theory on S 3 at strong 't Hooft coupling []19| . At 



zero and weak 't Hooft coupling, which is dual to a small AdS, thermodynamics of M = 4 
SYM theory on S 3 has been discussed in [p0j21 ]. In the free theory limit the Hagedorn 



and Hawking-Page temperatures coincide. At weak coupling it is not yet clear whether the 



transition is of first or second order [21]. Other studies of (Hagedorn) phase transitions in 



Yang-Mills theories include [11,12,22-35 



In this section we show that the critical Hagedorn behavior found in the last section 
arises also for a wide class of matrix quantum mechanical systems including M = 4 SYM 
on S 3 . Our discussion applies regardless of the order of the transition and also to strong 
coupling. In particular, we show explicitly that the Hagedorn divergences can be attributed 
to an effective potential of the form ( |2.31| ), which was only argued in the last section. 

The plan of this section is as follows. In next subsection we introduce the family of 
theories to which our discussion applies, which includes M = 4 SYM on S 3 . In the subse- 
quent subsections we discuss the large N expansion of these theories at finite temperature 
and identify new ingredients. We show that the free energy contain contributions from 
"vortices" and introduce a set of vortex diagrams to describe them. The vortex diagrams 
can be identified with degenerate worldsheets on the string theory side. The critical be- 
havior near the Hagedorn temperature and the effective action for the thermal scalar are 
recovered at the end. 
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3.1. Theories of interest 

Consider the following class of matrix quantum mechanical systems 

S= fair iVtr ^ (^(D T M a ) 2 - ^UqM^ + Mr^^(D T + £> a )£ a + V(M a , £ a ; A) 

(3.1) 

where: 

1. We have written the action in Euclidean signature, with the Euclidean time r having 
a period j3 = 4. In the zero temperature limit, /3 — > oo. 

2. M a and £ a are N x N bosonic and fermionic matrices respectively, and 

D T M a = d T M a - i[A, M a ], D T i a = d T t a - i[A, Co] • (3.2) 

are covariant derivatives. As a result, ( |Q| ) has a U(N) gauge symmetry, with M Q ,^ a 
transforming in the adjoint representation. The (0 + l)-dimensional "gauge field" A(r) 
plays the role of the Lagrange multiplier which imposes that physical states are singlet 
of U(N). M a ,£ a satisfy periodic and anti-periodic boundary conditions respectively 

M Q (r + /3) = M Q (r), £ a (r + /3) = -£ a (r) . (3.3) 



3. The frequencies w a and uj a in ( |3.1| ) are nonzero for any a and a, i.e. the theory has a 
mass gap and a unique vacuum. The number of matrices is greater than one and can 
be infinite. 

4. V(M a ,^ a ; A) can be written as a sum of single-trace operators and is controlled by a 
coupling constant A, which remains fixed in the large iV limit. 

M = 4 SYM on S 3 is an example of such systems with an infinite number of matrices 
when the Yang-Mills and matter fields are expanded in terms of spherical harmonics on 
S 3 . V(M a ,^ a ; A) can be schematically written as§ 

V = N (VXV 3 (M a , £ a ) + AF 4 (M Q , £ a )) (3.4) 



8 The precise form of the interactions depends on the choice of gauge. It is convenient to 
choose Coulomb gauge V • A = 0, in which the longitudinal component of the gauge field is set 
to zero. In this gauge, M a include also non-propagating modes coming from harmonic modes of 
ghosts and the zero component of the gauge field. 
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where V3 and V4 contain infinite sums of single-trace operators which are cubic and quartic 
in M a , £ a . A = gyM-N is the 't Hooft coupling. 

To end this subsection, let us recall the standard relation between the large N expan- 
sion of a matrix quantum mechanics like (|3.1|) (or a gauge field theory) at zero temperature 



with the string theory perturbative expansion In the large iV limit, the free energy 



of (|3.1|) can be organized in terms of the topology of Feynman diagrams 



00 

\ogZ = Y,N 2{1 - h) f h (\) (3.5) 

h=0 

where /o(A) is the sum of connected planar Feynman diagrams, and /i(A) is the sum of 
connected non-planar diagrams which can be put on a torus, and so on. The expansion 
( P . 5|) resembles the perturbative expansion of a string theory, with 1/N identified with 
the closed string coupling g s and //j(A) identified with contributions from worldsheets of 
genus-/i. For Af = 4 SYM theory on S 3 , //i(A) is the contribution of string worldsheets of 
h handles propagating in AdS$ x S5. 

In the next few subsections, we discuss the large N expansion of (|3.1| ) at finite tem- 
perature, and new ingredients arise. We find new contributions associated with Feynman 
diagrams with vortices, which can be identified with degenerate limits of a string world- 
sheet. As a result, the same critical Hagedorn behavior is recovered from gauge theories. 

3.2. Correlation functions in free theory 

In this subsection we discuss finite temperature correlation functions of (|3.1| ) in the 
free theory limit (i.e. with V = 0), focusing on the large N counting. We will find that 
at finite temperature, in addition to the standard 1/N 2 corrections due to non-planar 
diagrams, there are corrections due to vortices. This subsection makes preparation for the 
discussion of the interacting theory free energy in the next subsection. 

( [TTD has a U (N) gauge symmetry, which can be used to set the gauge field A(t) to 
zero. The gauge transformation, however, modifies the boundary conditions from (|3.3|) to 



M a {T + P) = UM a rf, Ca(r + p) = -Uti a tf . (3.6) 

The unitary matrix U can be understood as the Wilson line of A wound around the r 
direction (Polyakov loop), which cannot be gauged away. Correlation functions can then 
be written in terms of a path integral as 

(• " •>o,/9 = / dU j DM a {r)DUr) • • • e"*^-^ (3.7) 

15 



with M a , £ a satisfying the boundary conditions (|3.6|) . So and Z are the action and parti- 
tion function for the free theory respectively. The free theory action Sq has only quadratic 
dependence on M a and £ , thus the functional integrals over these variables in ( |3.7| ) can 
be carried out explicitly and (|3.7| ) can be reduced to a matrix integral over U only 



The free theory partition function can be written as [20,21 

Zo(P) = JdU e 1 ^ 

with Iq(U) given by 



Io(U) = J2-V n (0)TrU n TrU- 

71=1 



and 



z b (n[3) + (-l) n+1 z f (n[3), z b (J3) = J> 



-/3uj a 



(3.8) 



(3.9) 



(3.10) 



When the temperature T is small, the matrix integral fl3.8| ) can be evaluated in the large 
N limit as |2g£T | 



(3.11) 



where C is an ^-independent constant factor. Zq({3) becomes divergent if some V n {f3i) 
are equal to 1. From ( |3.1U| ) one can check that V\{(3) > V n (f3) for n > 1 and that V\{j3) 
is a monotonically increasing function of T, with V\((3 = oo) = and Vi(f3 = 0) > 1. 
Thus as one increases T from zero, there exists a Tjj, at which V\{Th) = 1 and Zq 
becomes divergent. Equation ( |3.11| ) only applies to T < Th- As pointed out in pp| , pT| , 
the divergence is precisely of the Hagedorn-type (|2.4p for a string theory in a spacetime 



whose Laplacian has a gap. The critical behavior of higher order terms in (|3.11|) near Th 
and the smoothing of the Hagedorn divergence at finite A" (i.e. in quantum string theory) 



for free Yang-Mills theory was further discussed in [11 



Correlation functions of gauge invariant operators can be obtained by first perform- 
ing Wick contractions and then evaluating the matrix integral for U. With boundary 
conditions (|3.6|), the contractions of M a and £ a are [[37 

$ab 

N 

(3.12) 



MS(r)M&(0) = -£ G s (T-mfcu a )U^ m U% 



m- 
oo 



&(r)&(0) 



Oab 

N 



{-l) m G f {T-m^a)U-rUT 3 



where G s and Gf are standard (0 + l)-dimensional propagators at zero temperature 

G s {r;u) = ±-e~ u H, G/(r;a;) = (-d T + ^)G s (r; u) . (3.13) 
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Fig. 4: An example of a double-line diagram at finite temperature. Each propa- 
gator carries a winding number (or image number), which should be summed over. 
Due to the presence of U- factors in ( 3.12|) , associated with each face one finds a 
factor of trU SA , instead of a factor N as is the case at zero temperature. 



It follows from (|3.12|) that at finite temperature, each propagator in a double-line Feyn- 
man diagram carries a winding number (or image number) m, which should be summed 
over (see fig. P . More explicitly, using (|3.12j ) , the contribution of a generic Feynman dia- 
gram to a correlation function of single-trace operatorsi can be written in the form0 |[7| 



/ 



N 2h- 



nn e 

i<j p =i m (p)__ 



\ 



V 



,<j P= i 




—trU Sl —trU S2 ---—trU 

N N N 



S F 



U 



(3.14) 

where i,j enumerate the vertices (i.e. operator insertions) and p enumerates the propa- 
gators between vertices i and j with the total number of propagators between them, 
label the images of G^ p \Tij). h is the genus of the diagram. In (|3.14|) , 



(...) 



u 



dU 



(3.15) 



with Iq(U) given by (|3.9| ). The powers s%, sz, • ■ ■ in the last factor of ( |3.14|) can be found 
as follows. To each propagator in the diagram we assign a direction and an orientation 
can be chosen for each face. For each face A in the diagram, we have a factor trlI SA , with 
s A given by 

s A = Y,(±)m ( *\ A=l,2,..-F (3.16) 



dA 



9 We assume the operators are normalized as Ntr( - • •). 

10 For notational simplicity, we only include bosonic modes in the equation below. It can be 
easily generalized to include fermions. 
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where the sum dA is over the propagators bounding the face A and F denotes the total 
number of faces in a diagram (see e.g. fig. In ( |3.16| ) the plus (minus) sign is taken if 
the direction of the corresponding propagator is the same as (opposite to) that of the face. 
sa has a precise mathematical meaning: it is the number of times that the Euclidean time 
circle is wrapped around by the propagators bounding a face A. We will thus call sa the 
vortex number for face A. Note that since the exponents of U add up to zero in ( |3.12|) , for 
each connected part of a Feynman diagram the sum of sa adds to zero. To illustrate more 
explicitly how (|3.14|) works, we give some examples in Appendix A. 

The partition function ( |3.8| ) and more generally matrix integrals in ( |3. 14j) can be 
evaluated to all orders in a 1/N 2 expansion. In Appendix B we prove that, up to corrections 
non-perturbative in N, the matrix integrals can be evaluated by treating each TvU n as an 
independent integration variable. More explicitly, ( |3.15| ) can be evaluated by replacing 



N 



1 



(3.17) 



i.e. 



where 



—TrU Sl — TrU S2 ■ ■ ■ — TrU SF 
N N N 



u 



From flgTTJ ), 



1 



n 

\i=l 



exp [-N 2 Y, v nW 



+ nonperturbative in N 



Vn(f3) 



1 " V n (0) 



— TrU Sl —TtU S2 
N N 



n=l 



n 




— TtU SF 


\ 


N 




( > 







+ 0(N 4 ) + nonperturbative in 



(3.18) 



(3.19) 



(3.20) 



where order 1/N 2 terms are obtained by contractions of one pair of </> Si 's, order 1/iV 4 
terms are obtained by contracting two pairs of 0's, and so forth. Each contraction brings a 
factor of jV 2„ 1 ■ Perturbative corrections in 1/N 2 terminate at order 1/N F (or l/iV F-1 ) 
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for F even (odd). For example, there is no other perturbative correction in 1/N 2 for the 
partition function (|3.11| ), and for F = 2 



1 

AT 



— trW—tvU" 1 ) = Sn^Smfi + T^5m+n,o + nonperturbative corrections . 



V\n\(P) 



(3.21) 



To summarize, combining (|3.14 ) and ( |3.2U| ) we find that for a correlation function of 
gauge invariant operators, there are two sources of 1/N 2 corrections: 

1. From the genus of the diagram as indicated by the power of 1/N in ( |3.14| ). This 
follows from the standard large iV counting. 

2. From the 1/N 2 corrections of the matrix integral ( |3.20|) . The leading order term in 
( |3.20|) imposes the constraint that for any face A of the diagram the vortex number 
sa should be zero. The next order corresponds to having nonzero vortex numbers in 
two of the faces, say the faces A and B with saSb 7^ and sa + sb = 0. Below, 
we will refer to those diagrams with nonzero vortex numbers as containing vortices, 
in anticipation of their interpretation from the string worldsheet0. From remarks 
below ( |3.16| ), if a face A of a Feynman diagram contains a vortex with vortex number 
sa, then the propagators bounding the face wrap around the Euclidean time circle 
sa times. At finite temperature, due to the presence of vortices, planar diagrams also 
contain higher order 1/A^ 2 corrections. 





Fig. 5: Examples of double-line diagrams with nonzero vortices. Each thin 
line (vortex propagator) represents a contraction in ( 3.20| ). Compare the left di- 
agram to fig. § Dia grams which are disconnected at zero temperature can be 
connected through vortex propagators as in the right diagram. 
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See also the discussion of [38] in the context of c = 1 matrix models. 
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It will be convenient to represents the vortex contributions diagrammatically: we represent 
each contraction in ( |3.20| ) by an oriented line between two surfaces which have the opposite 
vortex numbers. The orientation of a line is that it exists (enters) the surface if its vortex 
number is positive (negative). We associate a factor 1/N for each vortex and a factor 
l/v n (f3) to a line (vortex propagator) connecting two surfaces with vortex number ±n. See 
fig. |5] for some examples of such diagrams. Note that a diagram with otherwise disconnected 
parts connected by vortex lines should be considered as connected, as in the right diagram 
of fig. |^. In computing a correlation function one should sum over all possible vortex 
contributions. 

To summarize this subsection, in computing correlation functions at finite tempera- 
ture, one should consider not only Feynman diagrams which appear at zero temperature, 
but also diagrams with nonzero vortices. Explicit examples are given in Appendix A. 

3.3. Free energy in interacting theory and vortex diagrams 

We now consider the Euclidean partition function of the interacting theory. Our 
purpose is to identify Tjj and the critical behavior near Tjj to all orders in the 1/N 2 
expansion. 

In perturbation theory, the partition function can be evaluated by expanding the 
interaction terms in the exponent of the path integral 

00 (_l)n rP n 

z(p,\) = z (p)J2—r HdTityW'-vM^fi ( 3 - 22 ) 

n=0 U - J ° 1=1 

In ( |3.22| ), (• • -) a denotes free theory correlation functions and recall that V is given by a 
sum of single trace operators of the form Ntr( - ■ •). The free energy can be obtained from 

\ogZ = logZo + ~J~ / II dr < ■ ■■ V ^)) , p , connected (3.23) 

n=0 U - J ° i=l 

i.e. one sums only over the connected diagrams. The discussion in the last subsection for 
free theory correlation functions can now be directly carried over to logZ. In particular, 
there are two sources of 1/N 2 corrections: from the non-planar structure and from vortices. 
We can expand logZ in 1/N 2 as 

00 . 
log Z(J3) = N 2 - 2n Z n ((3) = N 2 Z ((3) + Z 1 ((3) + W MP) + ' • • (3-24) 

n=0 
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where Zq corresponds to the sum over connected planar diagrams with no vortices, while 
Z\ contains the sum of connected genus-1 non-planar diagrams with no vortices and planar 
diagrams with one pair of vortices, and so on. Recall that each vortex carries a factor 1/N 
and they always come in pairs. Also as remarked at the end of the last subsection, a 
diagram with otherwise disconnected parts connected by vortex propagators is connected. 




n 



Fig. 6: The propagators and vertices for vortex diagrams. The vertices ' n ' of a 
vortex diagram have n legs, each of which is labelled by a vortex number. The sign 
of the vortex number is positive (negative) if the corresponding leg exists (enters) 
the vertex. The total vortex number of a vertex is zero. We show Q(°' 2 \ Q^ 1 ' 3 ^ in 
the figure as illustrations. 



To elucidate the structure of Z g , we introduce a new set of "vortex diagrams", by 
generalizing the diagrammatical rules introduced below fig. [5]: 

1. Denote Q^ h,n ^ as the sum of connected Feynman diagrams with genus h and with n 
vortices. In terms of large N counting, Q < - h,n ^ is of order _/v 2_2/l_n , as we associate a 
factor 1/N with each vortex. Each vortex is labeled by a vortex number and the total 
vortex number carried by Q( h,n ^ is zeroS Diagrammatically, Q( h > n ^ are represented 
as vertices with n oriented legs. The leg exits the vertex if the corresponding vortex 
number is positive. 

2. Vortex diagrams are then constructed following the usual rules with Q^ h ' n > as funda- 
mental vertices and l/vb(j3),b > as propagators. Note that b is the vortex number 
carried by a propagator and Vf, was defined in (|3.19[) . 

3. The combinatoric rules are the same as standard Feynman diagram. In particular, 
if there are m identical vertices Q( h ' n ^ in a diagram, there is a factor 1/m!, which 
comes from the fact that disconnected diagrams are obtained from connected ones by 
exponentiation. 

12 This follows from the discussion below ( |3.f 6| ). 

21 



Using the above diagrammatical rules, we now enumerate the contributions to Z g . See 
fig. |6| for illustrations of propagators and vertices for vortex diagrams. 

Let us first look at Zq 1 which is given by the sum of all planar diagrams without vortex. 



In section 4 of |$7] it was shown that Zq is identical to the corresponding expression at 
zero temperature and thus is temperature-independentlll. Since the free energy — j3F is 
defined by subtracting the zero-temperature contribution (which is the vacuum energy) 
from ( p.24|) , we conclude that the planar contribution to the free energy is identically 
zerc@. 




Fig. 7: Vortex diagrams contributing to Z[ 



(3) 



We now look at Z\, which contains three contributions: (i) genus-1 contribution in 
free theory coming from the first term in ( |3.23| ); (ii) sum of genus-1 Feynman diagrams 
with no vortices; (iii) planar diagrams with vortices. The first contribution Zy is given 
by the logarithm of ( |3.11| ). The second contribution Z\ is given by QW. To find the 

(3) (0 2) 

third contribution Z\ , let us denote Q b ' the sum of all planar connected diagrams with 
two vortices of winding ±6. Graphically, it can be represented by a sphere with an arrow 
pointing in and an arrow pointing out, each carrying vortex number 6, as in the second 

(0 2*1 (3*1 

diagram of fig. |S|. Using Q b , Z\ is obtained by summing the vortex diagrams in fig. [7[ 



13 Zq is a special case of the discussion in section 4 of [37| with no external operator insertions. 

14 as is the case for a string theory below the Hagedorn temperature. 
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The combinatoric factor for a diagram with m vertices is 1/m following from the cyclic 
symmetry and we find that 



oo oo 



z?> = E E 



1 /«r 2) (A,/3) 



77?. 



6=1 m=l \ v ' / 6=1 

Adding all three contributions together we find that 



VbW) 



z x = + z[ 2) + z[ 3) 



= Q 1 >°(/?,A)-^(log(l 

6=1 

oo 



Ql°' 2) (A,/3 ) 



+ logv b ((3) 



Qi.0(/3, A) - $>g (v b {(3) - Qi°' 2) (A,/3) 



6=1 



(3.25) 



(3.26) 



Fig. 8: The dark thick line represents the resumed propagator Qh = 1 (0 2) 

v b -Q\ ' 



(3) (0 2) 

It should be clear from the above discussion of Z\ that Q h should not really be 
treated as a vertex. Rather all Q^' 2 ^ should be resumed along with the propagators -^rm 



to obtain a "resumed propagator" for each vortex number 



oo 1 

= E :=7* « 



(0,2)> 



)(0,2) 



(3.27) 



n=1^'"" ' V b -Q b 

as shown diagrammatically in fig. ^. Note that ( |3.26|) can be rewritten in terms of Q b as 



Z 1 = Q 1 '°(/3,A)+^log^(/3) 



(3.28) 



6=1 




Q 



(1.2) 



Q 



(0.4) 



+ 



Q 



(0.3) 




■ m - n n (0,3) 



Fig. 9: Vortex diagrams contributing to Z2. Compare plots in fig. || with the 
degenerate limits of a genus-2 surface in fig. Note the 2nd, 3rd and 5th diagrams 
in fig. H do not appear in the above since they contain propagators which have to 
have zero windings. 
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In the vortex diagrams for Z g with g > 2, only resumed propagators Qj, appear. As 
an example, the vortex diagrams contributing to Zi are shown in fig. |9|. Higher order 
diagrams contributing to general Z g can be similarly constructed. 

By now readers may have recognized the resemblance of vortex diagrams with the 
diagrams in fig. [I] and fig. [I| Indeed it is natural to identify vortex diagram contributions 
in the gauge theory with contributions from degenerate limits of string worldsheets in the 
corresponding string theory. For example, diagrams in fig. [)] can be identified with various 
degenerate limits (fig. ^ of genus two Riemann surfaces. In particular, vortices in gauge 
theory vortex diagrams can be identified with insertions of winding tachyon modes in the 
worldsheet. On the worldsheet if one follows a closed contour around the vertex operator 
of a winding tachyon mode of winding number b, the Euclidean time circle is traversed 
b times. Similarly, as discussed earlier if a face of a Feynman diagram contains a vortex 
with vortex number 6, the propagators bounding the face wrap around the Euclidean time 
circle b times. 

A more careful comparison between vortex diagrams for Z g and degenerate limits of 
a genus-g surface (e.g. between fig. |^ and fig. ||) also show some important differences: 

1. Notice that the 2nd, 3rd and 5th diagrams in fig. ^| do not appear in fig. |^. These dia- 
grams are distinguished in that some propagators are forced to have zero winding due 
to winding number conservation. One can convince oneself that this feature persists to 
all orders. Thus YM vortex diagrams do not correspond to the full contributions from 
degenerate limits of a Riemann surface. All propagators in the YM vortex diagrams 
carry nonzero windings. 

2. Various degenerate limits of a Riemann surface do not follow the standard Feynman 
rules and cannot be treated as Feynman diagrams. For example, the third diagram 
of fig. ^| can be obtained as a degenerate limit of the first diagram and the fifth as a 
limit of the fourth, etc. In contrast, the vortex diagrams we constructed in Yang-Mills 
theory do follow standard Feynman rules. In particular, different diagrams in fig. |9| 
do not overlap. 

Thus vortex diagrams correspond to a specific decomposition of the boundary of the moduli 
space and can be considered as defining an effective string field theory for the winding 
tachyon modes. 
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3. 4- Critical Hagedorn behavior and the effective action 

Now let us examine the critical Hagedorn behavior of (|3.24|) by increasing the tem- 
perature from zero. 

In free theory, as reviewed after equation ( |3.11| ), there is a Hagedorn temperature 
given by equation Vi(Ph) = 1 at which the free energy diverges as logZo ~ — log(/3 — (3h)- 
Note that there is only a one- loop divergence since all perturbative corrections in 1/N to 
( p. 11|) vanish. 

In the interacting theory the effective vertices Q^ h ' n > should be regular at any temper- 
ature since they involve only sums of products of ( ^.13| ) and their images. The divergences 
of Z n then can only occur when the resumed propagator Gb{0) ( |3.27|) become divergent, 
which happens when 

v b ((3) = Q ( ?> 2 \\,(3), i.e. =qQ> 2 \\,I3), 6=1,2-... (3.29) 

If we again assume that ( |3.29| ) is first satisfied for b = 1 as one decreases f3 from infinityEl, 
the Hagedorn temperature in the interacting theory is determined by 

Vi((3 H (X)) = 1 - Qi°' 2) (A, 0h(X)) (3.30) 

with the most divergent term in Z\ given by (see (|3.26|) ) 

Z 1 « - log(/3 - p H (\)) + finite, /3 ~ f3 H {\) . (3.31) 

Divergences in Z n can be analyzed following exactly the same power counting argument 
of the last section (after equation (|2.9|) ). We find that the most divergent contribution to 

Z n as (3 — > (3h is given by 

W^kr ■ (3 ' 32) 

Furthermore, since the construction of vortex diagrams follows the standard combinatoric 
rules of Feynman diagrams, we find that the most divergent pieces at each 1/N 2h order is 
precisely given by ([2.22 ) with the identification 



m 



I = v lW - Q r\x, 0), ^ = Qf^,-, + ' nM Qt%Q ( -i%a (3-33) 

15 which should be the case for A small since Q[ ' 2 ^ starts at order O(A). For large A, in principle 
this does not appear to be guaranteed from the gauge theory point of view. However, from string 
theory it appears always to be the case that the lowest winding modes become massless first. 
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where the subscripts in Q^ h ^ denote the vortex numbers for each leg. Similarly, one 
can use the same argument before ( |2.31|) to extract higher order terms in the effective 
action ( |2.31|) . The fact that we get (|2.31|) from divergences is guaranteed since vortex 
diagrams follow the rules of Feynman diagrams. It is also straightforward to work out the 
counterparts of ( |3.33| ) between A2/C in string theory side and Q^ m,n > . 

We note that on general grounds one expects that the free energy of the interacting 
theory can be written in terms of a matrix integral for U [ 2~I|] 

Z(p, A) = J dU e I(u) (3.34) 



with I(U) expanded in terms of all possible powers of trU r 



I(U) = Q (0) + Q ( n hrU n trU- n + £ Q^tvU^vU^vU 



77,^0 nmZ^O 

n-\-m-\-l = 

( 4 ) i.-Tm-i.^.TT'm+^.TTL 



+ J2 QnLp trUntrUmtrUltrUP + 



(3.35) 



n,m,l ,p^=0 
n-\-in-\-l -\-p=0 



where each Q[ n ) = YlhLo is a sum of contributions of diagrams of different genus 

h and • ■ • denotes windings of insertions. The vortex diagrams introduced earlier can be 
considered as the diagrammatical rules for computing ( |3.34| ). The effective action ( |2.31| ) 
then extracts the most important contribution near the Hagedorn temperature. 

It is important to emphasize that our discussion above should also apply to strong 
coupling. Q( n,m )(A), which are the basic building blocks of the vortex diagrams, can 
be defined non-perturbatively as follows. Since at each genus the number of Feynman 
diagrams grows with loops only as a power, we expect that Q( n,rn \\) should have a finite 
radius of convergence in the complex A plane. Once one obtains Q^ n,m \X) near the origin, 
one can analytically continue them to strong coupling. 



4. Conclusions and discussions 

In this paper we extracted Hagedorn divergences to all string loop orders and showed 
that they can be resumed. The resumed amplitudes have the form of an integral over the 
potential ( |1 . 1[ ) for the thermal scalar and smooth the divergences. We presented arguments 
both from a worldsheet approach and from Yang-Mills theories using AdS/CFT. In the 
double scaling limits (|2.26| ), worldsheets with arbitrary number of thermal scalar insertions 
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become equally important, which is consistent with the expectation that the thermal scalar 
will condense and the spacetime background will shift. 

The fact that one can obtain the thermal scalar potential to arbitrary higher orders 
by analyzing the local divergences in the thermal string phase is interesting. The potential 
would enable one to find other possible phases of the theory. The results also give an 
unambiguous prescription for computing the potential for the thermal scalar near the 
Hagedorn temperature from string amplitudes. The relation we found between vortex 
diagrams in Yang-Mills theory at finite temperature and degenerate limits of worldsheet 
Riemann surfaces is rather intriguing and worth investigating further. 

Finally we note our strategy for extracting the thermal scalar potential should also 
be applicable to the tachyon condensation in a circle with anti-periodic boundary condi- 
tions (for a recent discussion see p9[]). 
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In this appendix we give some explicit examples on the use of equation ( |3.14j) for cal- 
culating correlation functions between single trace operators. For definiteness we consider 
only bosonic operators, but the procedure is analogous for operators involving fermions. 
Consider the following simple example 



(NTtM 4 (r)NTrM 4 (0)) 



(A.l) 



where M can be any of the bosonic modes in ( |3.1| ). The calculation of ( |A.1|) amounts 
to drawing all possible double line diagrams. For example the disconnected planar con- 



tribution is given in fig. [L0|. From (|3.12| ), each propagator carries an image number (or 
winding number), which should be summed over. Each face A carries a factor trlI SA . sa 
is determined by choosing a direction for the propagators, and an orientation for the face, 
as explained below ( |3.16| ) . fig. [10| therefore gives a contribution of the form 



^2 J2 G s (~mP)G s (-nP)G s (-pP)G s (-qP)(TrU m TrU n TrU- m - n TrU p TrU q TrU- p - q ) 

m,n,p,q= — oo 

(A.2) 



u 



j-j—m-p-n+q 




TiU q - m * 




r-m-p+q 



Fig. 11: Planar connected contributions to (TrM 4 (r)TrM 4 (0)) 




Fig. 12: Some non-planar (torus) connected contributions to (TrM 4 (r)TVM 4 
For visualization purpose, the edge of one of the faces is drawn in red. 
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The connected planar contributions are given in fig. 11 with, for example, the first 
diagram given by 



G s (t- m[3)G s (T - n[3)G s (r - p(3)G s (r - q/3) (TrU m - n TrU n - p TrU p - q TrU q - rn ) 



N 2 

m,n,p,q= — oo 



U 



(A.3) 

In fig. |12] we have also plotted some connected non-planar diagrams, with the first diagram 
given by 



u 



G s (r-m/3)G s (T-n/3)G s (T-p/3)G s (r-q/3)(TrU m -P +q - n TrU- m+ P- q+n ) l 

m,n,p,q= — oo 

(A.4) 




Fig. 13: Connected vortex diagram from disconnected double line diagram 



Now let us consider the evaluation of the expectation values of traces of U in ( |A.2j )- 
( |A.4p using ( pj.2U| ). At leading order in the large N expansion the expectation values give 
N F , where F is the number of traces, times some product of Kronecker delta enforcing all 



exponents to be zero. In this case we recover the results of [p7| . Higher order corrections 
in 1/N can be described graphically by inserting pairs of vortices on different faces of the 



diagrams and connecting them with the propagator 



MP)' 



One should sum over all the 



possible ways of inserting pairs of vortices. Note that each vortex insertion gives a factor 
of 1/N. Diagrams with disconnected parts connected by vortex propagators should be 



considered as connected as in fig. 13. Note that in terms of large iV counting fig. 13 is of 



the same order as those in fig. 12 with no vortices 
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Appendix B. Proof of ( gTT8|) 



In this appendix we prove equation ( |3.18| ). In the next subsection we discuss some 



elementary aspects of U(N) group integrals. We then proceed to evaluate (|3.8|). Equation 
( p.!8|) is proved in the end. 

B.l. Group integrals over U(N) 

Consider the following integral over the unitary group U(N) 

I = ^- j dU \\{TiU ai ) bl f[(TrU-^ (B.l) 

i=l j=l 

where Oj, bi, q, di are positive integers and 

k s 

D = oA = c i d J ■ ( B - 2 ) 
i=l j=i 

Vn is the volume of U(N). 

Products of traces of U can be expanded in terms of characters of irreducible represen- 
tations of U(N), which are in one to one correspondence with irreducible representations 
of the symmetric group (see for example [f40|), 



JJiTrU^h = J2xx(ai,bi)xx(U) (B.3) 



i=i 



where A labels the irreducible representations of the symmetric group So- X\( a i>bi) is the 
character of the conjugacy class0 of So given by the set {(a^, bi)} in the representation A. 
Xx{U) is the character of U in the irreducible representation of U(N) labelled by A. Now 
by using the orthogonality property for characters we can write: 

I = Y l Xx(nA)xA*,<li)^- f dUxx(U)xx'(Jf) 

VN 1 (B.4) 
= ^2xx(ai,bi)x\(ci,di) . 



16 Recall that two elements of Sd are conjugate if and only if they consist of the same number 
of disjoint cycles of the same lengths. Denote the number of cycles of length en by bi then a 
conjugacy class in Sd is given by a set of k couples {(en, bi)} i = 1, ..k such that ^ fe =1 ciibi = D. 
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The evaluation of (|B.4|) can be divided into the following two cases: 



1. If D < TV, then the sum over A can be evaluated giving [40 



I = <*{(oi,&0}.{(ci.di)} Yl Xx ^ hi ^ 2 = S {^iA)},{(a,d i )} Y[^ bl W (B.5) 

A i=l 

where the completeness of characters of the symmetric group Sp enforces the sets 
{(aj, bi)} and {(q, dj)} to define the same conjugacy class in Sd, i-e., to be the same 
apart from reordering. This means that the integral is zero for D < N unless for any 
factor of Tr[U a ] b in the integrand there is a corresponding factor of Tr[U~ a ] b . 
2. If D > N one needs to restrict the sum over the irreducible representations A to the 
representations where X\(U) ^ 0, that we will indicate formally as A < N. In this 
case the result is more complicated and we do not have a closed form expression. For 
the case in which the sets {(a^, bi)} and {(c^, di)} are equal up to reordering one has 

k 

1= J2xx(ai,bi) 2 <]Ja i b %\ . (B.6) 

\<N i=l 

B.2. Partition function integrals 

We now consider the evaluation of the free theory partition function (|3.8|) . To warm 
up let us consider the following integral 

/ dUe z ^ UTrU ' = -i- / dUY i(^TrC/TrC/t)P 
Vn J V N J ^ pi 

N 

= Y,4 + o(z?) ( B - 7 ) 

p=0 

For < z\ < 1 the corrections to the N = oc result are of order O(z^) and are therefore 
exponentially suppressed in N. In the more general case ( |3~8l ) (with V n (f3) = z n ) one can 
proceed exactly as above, writing 



= _L [ dUe 10 ^ = -i- [dUexply ^TrCTTrC/t' 
V N J V N J * n 

\n=l 

1 f ^ / „2L 7 Pr,. \ 

= — / dU TT -4 1 — (Trt/ n Trt/" n ) Pn (B.8) 

n=l \p n =0 / 

oo 1 
= IIl^-^) 
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where C(N) is given by 



C(N) 



oo / oo \ 1 /" oo/cxd p 

n E*U/*n 

.n=l Vp„=0 / KArJ n =l \p n =0 P "- R 



(Tr£/ n Trt/- n )^ 



< 



oo / oo 

n e 

,n=l \p„=0 



r»(iv) 



V np n >N 

Z — ^ n 



V n Pn >N 

Z — 



(B.9) 



Note the the subscript in the above equation indicates that one should only sum over those 
p n which satisfy ^ n np n > N. D(N) can be estimated as follows. Consider the expansion 



oo ^ oo 

II 1 _ VTn = Z)°n(^>^»-)^ 



1 — Z n t 

n=l n n=0 

where a n are polynomials in the Zi with positive coefficients. Note that 



(B.fO) 



Define 



D(N)= E 

n=N+l 



z* — max(^i, z^ i z^ , Zn , ...) 



(B.11) 



(B.12) 



Below Th, we have z* < 1. Then we have < a n (zi, Z2, ■■■) < a n (z*, zj, z%, ...) = b n z™ 
where the 6 n 's are the coefficients of the series of IIm=i i-z m = S^Lo ^n-z* • This series 
has radius of convergence equal to f because the function has no singularities for \z*\ < 1. 
It then follows that for a given e > there exists an M(e) such that for n > M(e) it is 
true that b n < (I + e) n . Then for e < j- - 1 and N > M (e) the following holds: 



c(n) < d(n) < J2 ((i+^o n = ^r?r * J 

n=AT+l 



(1 + e)z* 



(B.13) 



and therefore the corrections are exponentially small since (1 + e)z* < 1. 
To summarize we find that 



7 Q = _L / dUexp ( Y ^TtU n TrU^) = f\ — 



-Nc 



(B.14) 



where c = - log(z*) > and K > 0. 
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B.3. Correlation functions 

Correlation functions (|3.15 

k 



Y^(TrU ai ) bi JJ(TrJ7 _c ') d ' \ = 4~ / dU e Iu{u) f[(TrU a *) h f[(TrU- c i) d > (B.15) 

i=l j=l / [/ ° *=1 5=1 

where cii,bi,Ci, di are positive integers of order O(N ) can now be calculated easily using the 
technique above. Correlation functions of the form (y[ n (TrU an TrU~ an ) bn ^ v are obtained 
by taking derivatives on Z ( |B.14[ ) with respect to ^ 



\ 1 r/^™ 7 
Y\(TrU an TrU- an ) K ) = — T] n bn ^ . (B.16) 

If in ( [B.15Q the {(aj,6j)} are not matched with {(cj,<ij)} up to reordering, due to (|B.5| ), 
the correlation function is zero up to nonperturbative corrections which are of order (z*) . 
For example, (TrU a TrU a TrU- 2a ) u is zero at any finite order in i expansion unless a is 
zero. 

The above results can be summarized by the following: the integrals can be evaluated 
by treating each TvU n as an independent integration variable. More explicitly, replacing 

llUT-^, ^W"^ <P-n = C, 0o = l (B.17) 



then 



— TrU Sl —TtU S2 ■ ■ ■ — TrU SF 
N N N 



u 

oo 



z. 



1 f / 1 — Z \ 

— / J^J <i</>i<i</>* <p Sl ■ ■ -(/} SF exp | —N 2 -0 n 0* J + nonperturbative in N . 

(B.18) 
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